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We present a search for anomalous semileptonic decays of heavy flavor hadrons produced in 
association with a W boson, in proton-antiproton collisions at -^8=1.96 TeV. We use 162 pb“^ of 
data collected with the CDF II detector at the Fermilab Tevatron Collider. We select events with 
one W boson and at least one jet with an identihed secondary vertex. In the jets with a secondary 
vertex we look for a semileptonic decay to a muon. We compare the number of jets with both 
a secondary vertex and a semileptonic decay, and the kinematic properties of these jets, with the 
standard model expectation of W plus heavy flavor production and decay. No discrepancy is seen 
between the observation and the expectation, and we set limits on the production cross section of a 
B-like hadron with an anomalously high semileptonic branching ratio. 

PACS numbers: 13.85 Qk, 13.38 Be, 13.20 He 


Anomalous semileptonic decay signatures are pre¬ 
dicted by some theories of physics beyond the standard 
model (SM) P|. A light supersymmetric partner of the 
bottom quark in a sample of ordinary b quarks could, 
for instance, enhance the apparent semileptonic branch¬ 
ing ratio of B decays. In this scenario, the semileptonic 
decay of supersymmetric B mesons, characterized by lep¬ 


tons and charmed mesons B^ivDX or B^iDX 
(where L indicates a scalar neutrino and X labels possible 
additional hadrons), could look similar to standard model 
B decays. In a sample of events with a W boson and at 
least one jet in Run I of the Tevatron Collider, CDF ob¬ 
served an excess of events in which a single jet contained 
both an identified secondary vertex and a muon • The 
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rate of these events was inconsistent with that expected 
based on the semileptonic branching fractions of heavy 
flavor jets (6 or c). In this analysis we follow up on this 
observation with the larger data sample of Run 11. The 
CLEO Collaboration has performed a direct search for 
a supersymmetric B meson with mass between 3.5 and 
4.5 GeV/c^, just below the e+e“ ^ BB threshold, find¬ 
ing no evidence for a light scalar bottom quark Q. The 
D0 Collaboration has also recently reported on a search 
for excess of semileptonic decay in events with a W bo¬ 
son _Q, setting upper limits on the production of non-SM 
Wbb-\ike and top-like processes. 

We present the analysis of a sample of events with as¬ 
sociated production of W boson and quarks or gluons 
(referred to as IT-I- jets in this article). We identify the 
W boson by its decay to an isolated electron(muon) car¬ 
rying large transverse energy (momentum) Et(j)t) with 
respect to the beamline, plus a neutrino. We refer to 
these high px electrons or muons as “primary leptons”. 
The neutrino escapes the detector and results in an im¬ 
balance of total transverse energy, referred to as “missing 
Et^^t) 0. Quarks produced in association with the 
W boson hadronize and are observed as jets of charged 
and neutral particles. We identify the heavy flavor (h.f.) 
jets from their displaced secondary vertices. Indepen¬ 
dently, we search for semileptonic decays of h.f. hadrons 
by looking for a low momentum muon within the jet. At 
the Tevatron, the primary SM contributions to W-l-h.f. 
quarks in the final state arise from the production of W66, 
Wcc, Wc and tt, with small contributions also from WZ 
(with Z —> 55, cc) and single top quark production. The 
dominant background is due to QCD production of mul¬ 
tijets, where one of the jets is misidentified as the primary 
lepton, and there is large.^T in the event. Moreover, one 
of the jets must contain a h.f. quark or be misidentified 
as containing one, referred in this article as “mistagged”. 
Smaller background sources arise from W plus a misiden¬ 
tified h.f. jet, Zbb, Zcc^ WW, ZZ and Z ^ tt events. 

The strategy of this analysis can be summarized as 
follows: We start from the W-l-jets data sample with at 
least one jet having a displaced secondary vertex. The 
same sample has been studied and used to measure the 
tt production cross section . In this sample we find the 
number of events in which at least one jet has both a dis¬ 
placed vertex and a muon, and compare it to expectation. 
The expectation is estimated from data and Monte Carlo 
(MC) simulation, determining the probability of finding 
a muon within the jet with a displaced vertex. This 
probability depends primarily on the semileptonic decay 
branching ratio (B.R.) of the R, the kinematics of SM 
processes that lead to W-l-h.f. event, and the efficiency 
of muon identification. The expectation determined with 
such an approach is relative to the number of observed 
events with a displaced secondary vertex, hence we effec¬ 
tively compare the observed rate of semileptonic decay 
to muons with that expected from the SM processes, and 
use this comparison to set a limit on R-like decays with 
high semileptonic decay B.R. In addition, sensitivity to 


new physics with SM-like semileptonic decay B.R., and 
hence limits on such processes, can be calculated from 
the comparison of the observed and predicted number of 
events whose jets have a displaced vertex only, without 
further requiring a semileptonic decay. 

The CDF II detector is described in detail else¬ 
where 0. The components used in this analysis in¬ 
clude the silicon vertex detector (SVX), the central outer 
tracker (COT), the central electromagnetic and hadronic 
calorimeter, and the central muon detectors (CMU, CMP 
and CMX). The data sample, produced in pp collisions 
at -^5=1.96 TeV during Run II of the Fermilab Tevatron, 
was collected by the CDF II detector between March 2002 
and August 2003. This analysis is based on an integrated 
luminosity of 162±10 pb“^ (150±9 pb“^ with the CMX 
detector operational MC simulations are performed us¬ 
ing PYTHIA v6.2 11 and HERWIG v6.4 |i[l3, cou¬ 
pled with CTEQ5L parton distribution functions m 
Mode ling of 5 and c hadron decay is provided by QQ 
v9.1 [ij. Diboson and W-l-h.f. production are instead 
modeled with ALPGEN which generates high mul¬ 
tiplicity partonic final states using leading order matrix 
elements, and passed to HERWIG for parton showering. 

A detailed description of the event selection is given in 
0. Differences specific to this analysis are due mainly 
to small improvements in the calibration of the muon 
chambers and the restriction of the primary muon accep¬ 
tance to the best understood regions of the detector, as 
done in |I4|. Candidate primary electrons(muons) from 
W decay must have Et{pt) greater than 20 GeV with 
|ry| < I.l, and the event must have.^T > 20 GeV, consis¬ 
tent with the presence of a neutrino from the W boson 
decay. Since leptons from W decay are expected to be 
isolated from other energy deposits, the isolation (/), de¬ 
fined as the calorimeter transverse energy in a cone of 
All = \J{Arjy + (A(/))2 = 0.4 around the lepton (not 
including the lepton energy itself) divided by the lepton 
transverse energy/momentum, is required to be less than 
0.1. Finally, the W-l-jets dataset is partitioned according 
to the number of jets with Et > 15 GeV and |? 7 | < 2.0 in 
the event. Jets are identified using a cone algorithm with 
a cone opening of 0.4 and are constrained to originate 
from the pp collision vertex. The jet energies are cor¬ 
rected to account for variations of the detector response 
in r] and time, as well as multiple pp interactions. 

Heavy quarks (5, c) are identified by the presence of 
secondary decay vertices reconstructed with silicon-based 
tracking information. The search for a secondary ver¬ 
tex is performed jet by jet. An event by event primary 
vertex nearest to the identified lepton from the W de¬ 
cay is first selected. After identifying a secondary ver¬ 
tex, the two-dimensional decay length L 2 D is calcu¬ 
lated. Secondary vertices corresponding to the decay of 
h.f. hadrons are expected to have large positive L 2 d- A 
jet with L 2 D greater than 3 times its uncertainty is said 
to have a “SecVtx tag”. The average efficiency for tag¬ 
ging 5-jets is a function of the jet Ex and ranges from 
about 20% for Ex — 20 GeV to around 40% for jets 
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with Et above 40 GeV 0. Table HI shows the number 

TABLE I: Total number of IT+jets events with at least 
one SecVtx tagged jet, and the estimated sample composi¬ 
tion according to the SM contributions. The rightmost col¬ 
umn shows the average efficiency of additionally tagging the 
SecVtx-tagged jet with the SLT. 


W+ n jets [n]: 

1 

2 

> 3 

eSLT [%] 

Total Data 

157 

70 

59 

- 

Wbb 

37T11 

19±6 

6.0T1.6 

7.1T0.7 

Wcc 

14±3 

6.9T1.9 

2.0±0.6 

4.1T0.5 

Wc 

34±9 

6.7T1.7 

1.5±0.4 

4.1T0.5 

Zbb 

0.7T0.2 

0.3T0.1 

0.09±0.09 

7.1T0.7 

Zee 

0.4±0.1 

0.3±0.1 

0.09±0.09 

4.1±0.5 

WW 

0.6T0.1 

1.0T0.2 

0.41T0.06 

3.1±0.5 

wz 

0.5T0.1 

1.0T0.2 

0.31T0.03 

6.2±1.3 

zz 

O.OiO.O 

O.ITO.O 

O.OOiO.OO 

8.8±2.7 

Z ^ TT 

1.1T0.2 

O.lTO.l 

O.OOiO.OO 

3±4 

QCD 

24±3 

9.2±1.7 

5.0±1.1 

11±2 

W+l.f. 

41±6 

15±2 

7.1±0.8 

1.2±0.2 

top (s-chann.) 

0.6±0.1 

1.8±0.2 

0.7±0.1 

7.8±0.8 

top (t-chann.) 

2.0±0.2 

2.2T0.3 

0.6±0.1 

7.0±0.8 

ti 

0.9±0.1 

6.1T0.5 

35±3 

8.4±0.8 


of events passing our selection requirements with at least 
one SecVtx-tagged jet. We also show here the expected 
composition of the observed SecVtx tagged sample as 
determined by the same techniques used in the ti cross 
section analysis As in that paper, the tt contribu¬ 
tion is taken to be the difference between the observed 
number of events and the total SM expectation, exclud¬ 
ing tt itself. The corresponding measured ti cross sec¬ 
tion is 5.3±1.4 pb. The small component from ti in the 
W-l-1,2 jets sample is calculated using the expected ratio 
of reconstructed ti events in W+ > 3 jets, to W+1 and 
W+2 jets. Finally, for W+1,2 jets, we scale the overall 
sample composition determined in [^, and expressed in 
terms of expected tagged events, in proportion to the ob¬ 
served number of tagged events in this analysis. Hence, 
the sum of the components in Table H] is equal to the 
number of observed events by construction. The scaling 
between expected and observed events differs from unity 
by approximately 10% or less. In Tabled QCD refers to 
non-top QCD multijet production, and VF-l-l.f. indicates 
W plus light flavor jets, where one or more jets has been 
wrongly tagged by the SecVtx algorithm. 

Muons in jets are identified by matching the tracks 
of the jet, as measured in the COT, with track seg¬ 
ments from at least one of the three muon detectors. 
The distance between the extrapolated track and the 
muon chamber track segment depends upon muon multi¬ 
ple scattering and, at high px, intrinsic resolution of the 
drift chambers. A single quantity (L), with the distribu¬ 
tion centered around zero and normalized by its expected 
width, is constructed using parametrization of the match¬ 
ing resolutions directly measured from J/ij} —> 

Z W and W ^ pLv data M- Candidate tracks 


are identified as muons if \L\ < 3.5, where the effi¬ 
ciency for identifying a muon plateaus at about 85%. 
We refer to such muons as “Soft Lepton Tags”or SLT. 
Backgrounds are preferentially rejected because they are 
mostly accidental overlaps of uncorrelated tracks and 
hadrons that reached the muon chambers, for which 
the track-stub match is notably worse than for muons. 
Muons with px > 5 GeV/c are selected and required 
to be within Ai? < 0.6 of a jet axis. The probability 
of misidentifying a hadron as an SLT muon, denoted as 
the “mistag” probability, is measured from data using a 
large sample of jets in events collected with an inclusive 
photon trigger (y-l-jets sample). The number of mistags 
is counted on a per-track basis for jets with Ex > 15 
GeV, and an SLT mistag probability look-up table is con¬ 
structed as a function of the track px, rj and (j). The av¬ 
erage mistag probability is approximately 0.7% per track 
for tracks above px = 5 GeV/c. The tags in the y-l-jets 
sample include 21.0 ± 1.4% [l^ that are due to genuine 
h.f. decay. We subtract this fraction from all mistag 
predictions based on the mistag probability look-up ta¬ 
ble. The number of tags predicted with this technique 
agrees to within 10% |l4l of the observation in a variety 
of samples, including QGD multijets and ^-l-jets events. 

We determine the expected number of events with at 
least one jet tagged by both SLT and SecVtx (referred 
to as “SLT-SecVtx tag”) combining the contributions 
in Table m with the efficiency (cslt) for finding an SLT 
muon tag in the same jet tagged by SecVtx. For a h.f. 
jet, esLT depends primarily on the semileptonic B.R., 
B ^ p v X, and the requirement of 3 GeV/c as the 
muon minimum px, but it is also affected by the over¬ 
all efficiency of SLT tagging and correlations with the 
selected jets. For instance, measured jet energy is under¬ 
estimated for jets containing muons and neutrinos, and 
consequently low Ex jets with semileptonic muon decays 
can fail the Ex > 15 GeV requirement. This effect is 
well modeled in the simulation. Since esLT depends on 
kinematic properties of the jets and muons, it is there¬ 
fore specific to jets from different physics processes. For 
Whi), Wcc, ti, single top, diboson events and Z tt, 
the kinematics are determined from MG simulation and 
the efficiency for given kinematics from the data. In sim¬ 
ulating the SLT tagging, we use both the efficiency of 
tagging a muon in the jets, measured from data, and 
the probability of tagging a jet by accidentally tagging a 
hadron, described by the mistag probability look-up ta¬ 
ble. The average values of cslt are shown in Table ||| As 
a cross check of this technique for predicting the num¬ 
ber of SLT-SecVtx tags in VF-l-jets events, we apply the 
same procedure to a high purity bb sample, derived from 
SecVtx-tagged dijets and whose composition is measured 
to be ~ 95% bi) and ~ 5% cc 0. The number of ob¬ 
served SLT-SecVtx tags in this sample agrees with the 
prediction to within 10%. 

The value of cslt for the QGD background component 
is determined directly from data. We select events with 
non-isolated (/ > 0.15) primary lepton, in order to pref- 
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erentially select QCD events, and require > 20 GeV. 
In QCD events, high often results from a neutrino 
from a semileptonic decay in a h.f. jet, or when a jet 
is pointing towards one of the non-instrumented regions 
of the calorimeter. The latter typically correspond also 
to gaps in the absorber material through which hadrons 
can escape, creating false signals in the muon chambers. 
Jets separated from the primary lepton are required to 
have a SecVtx tag, and we measure cslt in these jets to 
be 11±2%. This value is independent of the jet multi¬ 
plicity, the primary lepton’s flavor and isolation, while 
it increases with higher Finally, the probability of 
tagging with the SLT a light flavor jet that has been 
wrongly tagged by the SecVtx algorithm is determined 
using a sample of jets with negative L 2 D in VF-|-jets MC, 
using the track multiplicity of the jets and the SLT mistag 
probability. 

The number of VF-|-jets events with at least one SLT- 
SecVtx tag is shown in Table HIl together with our expec¬ 
tation from SM processes. There is no excess of observed 
events with SLT-SecVtx tags. The quoted uncertainty in 
the SM expectation is the total systematic uncertainty. 
Since we effectively measure the ratio of events with SLT- 
SecVtx tagged jets to all SecVtx-tagged jets, systematic 
uncertainties in the event selection, acceptance and re¬ 
construction efficiency cancel out. Systematic uncertain¬ 
ties considered here include the efficiency for SLT tag¬ 
ging muons inside jets from heavy quark fragmentation 
(-1-1/ — 8%), or from mistags (±10%). The uncertainties 
also include a contribution from MC modeling of particle 
tracking efficiency in dense jet environments (±5%), and 
a contribution from the MC modeling of the jet energy 
calibration (±2%), which enters as a correlation between 
the jet Et > 15 GeV requirement and the presence of 
a muon in the jet. Uncertainty on the SM semileptonic 
decay B.R. of B mesons (±4.4%) is also included. 

We compare the kinematic properties of the SLT- 
SecVtx tagged sample with SM expectations. Figure Q] 
shows the distribution of L 2 D (a) and Et (b) for the 
SLT-SecVtx tagged jets, and pr (c) and (d) for the 
SLT muons. The quantity is the component of the 
muon momentum orthogonal to the jet axis, which is a 
discriminant for h.f. decays where the muon can receive 
a large boost orthogonal to the jet axis. Kolmogorov- 
Smirnov tests of the shape of the distributions indicate 
good agreement with the SM predictions. Figures ^^c) 
and Ed) have one extra data entry because one jet has 
two identified SLT muons. The shape of the kinematic 
variable distributions from the SM processes are obtained 
from MC simulations, with the exception of the QCD 
component, which is derived from data with the primary 
lepton isolation I > 0.15. 

Since semileptonic decays in jets produced in associa¬ 
tion with a W boson are in agreement with the SM rate 
and kinematics, we proceed to set a limit on non-SM pro¬ 
cesses that may lead to an excess of SLT muons. We first 
derive a model-independent limit on the excess of events 
above SM, as a function of the jet multiplicity. A limit 


TABLE II: Number of expected and observed events with at 
least one jet tagged by both the SLT and the SecVtx algo¬ 
rithm. The total expected is based on the contributions from 
all SM processes. The 95% C.L. limit is calculated for the 
number of events above SM. 


W+ n jets [n]: 

1 

2 

>3 

Wbb 

2.63±0.80 

1.56±0.45 

0.47±0.13 

Wcc 

0.60±0.15 

0.27±0.08 

0.08±0.03 

Wc 

1.51±0.40 

0.26±0.07 

0.06±0.02 

Diboson, 




Z ^ TT 

0.09±0.06 

0.12±0.02 

0.03±0.01 

Z+bb, cc 

0.07±0.01 

0.03±0.01 

0.01±0.01 

W+U. 

0.70±0.13 

0.20±0.04 

0.08±0.01 

QCD 

2.70±0.62 

1.00±0.26 

0.55±0.16 

top (s+t channel) 

0.19±0.02 

0.33±0.03 

0.11±0.02 

tt 

0.11±0.01 

0.69±0.07 

3.56T0.43 

Total Expected 

8.6±1.2 

4.5±0.8 

4.9±0.4 

Observed 

7 

5 

1 

Excess 95% C.L. limit 

6.7 

6.8 

3.5 



(c) (d) 


FIG. 1: Kinematic distributions for SLT-SecVtx tagged jets. 


for specific models can be obtained by introducing the 
model-dependent SLT tagging and event reconstruction 
efficiency. We derive limits on two kinds of models, one 
with the signature and kinematics similar to that of SM 
Wbb (1F65-Iike), and one that is similar to top pair pro¬ 
duction (tt-like). In each, we introduce a generic model 
where the non-SM iJ-like object, B, has all the kine¬ 
matic properties of a B hadron when it decays semilep- 
tonically, but decays with the semileptonic B.R. higher 
than the SM value. The limits follow from the model- 
independent case by applying the different SLT tagging 
and event reconstruction efficiency for the two cases. We 
set 95% C.L. limits using a Bayesian approach, with the 
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likelihood given by a Poisson distribution, multiplied by 
a Gaussian distribution restricted to positive values de¬ 
scribing the background, and a constant probability den¬ 
sity function representing the signal cross section. Table 
im gives the model-independent 95% C.L. limits for the 
number of VP-|-jets events with a SLT-SecVtx tagged jet, 
coming from anomalous semileptonic decays. 

Limits on production cross-sections can be determined 
by considering the model-specihc cslt, the event recon¬ 
struction efficiency and the luminosity of the data sam¬ 
ple. Acceptance and efficiency are calculated for 1 and 
2 jet events together in the case of a Whh-\\ke process, 
and for >2 jets in the case of a tt-like process. As al¬ 
ready noted, the efficiency for hnding an SLT muon in a 
SecVtx tagged jet depends on the leptonically decaying 
species present, their branching ratios and the jet kine¬ 
matics. In semileptonic B meson decay, a fraction fd of 
detected muons originate from the direct decay of the b 
(with B.R. 10.68±0.22% M) , while the rest are from 
the sequential decay b ^ c ^ ^ v X. The total SLT 
tagging rate depends on the efficiency (e^) of tagging 
the primary and secondary decay muons, with an addi¬ 
tional contribution from mistags {cm)- For Wbb events 
we estimate = 7.1 ±0.7%, with = 5.9 ±0.5% and 
fd = 56±3%. For ti events we estimate = 8.3±0.8% 
with = 6.8 ± 0.6% and fd = 58±2%. For a non-SM 
particle B with high semileptonic B.R., the efficiency of 
tagging with the SLT a SecVtx-tagged jet can therefore 
also be expressed by: = k ■ fd ■ + Cm, where k 

parametrizes the unknown semileptonic B.R. of the R, 
and we assume no contribution from sequential decays. 
The possibility of multiple SecVtx tagged jets is properly 
taken into account in the calculation of the efficiency, 
e^LTi for an event with one or more SecVtx tagged jets, 
to have that same jet additionally SLT tagged. The up¬ 
per limit on the number of non-SM events present in the 
SecVtx tagged sample, before requiring SLT tags, is ob¬ 
tained by including in the Poisson distribution for 
the signal: 

P(s|n) = ^SLT + ^) ^ 

n\ 

where b is the background (i.e. the total SM contribu¬ 
tion), s is the non-SM signal before requiring SLT tags, 
and n is the number of observed SLT-SecVtx tagged 
events. Since the non-SM events s must be a subset of 
the observed events shown in Tabled the background in 
the likelihood must also decrease proportionally to the 
increase of s. We introduce the average SM efficiency 
esM = b/Ntot, where Wot is the total number of events 
with a SecVtx tag but before SLT tagging, so that the 
numerator of Equation ^ and the argument of the expo¬ 
nential are replaced by: 

^ + esM(Wot — s) = s ■ ~ ± b. (2) 

In Equation|21 the SM contribution decreases as a func¬ 
tion of s and the sensitivity of the measurement exists 



FIG. 2: Upper limits on the production cross section of Wb&- 
like and tt-like processes as a function of the branching ratio: 
B.R.(R ^ X). 

as long as the average efficiency of the non-SM process 
differs from the efficiency averaged over the SM contri¬ 
butions. This limitation is a direct consequence of the 
fact that this analysis measures the ratio of SLT-SecVtx 
tagged over all SecVtx tagged events. 

Finally, the event efficiency, including the geometric 
acceptance, the set of cuts applied to the primary lepton 
and to the jets, as well as the efficiency of the SecVtx 
algorithm to tag h.f. jets, is measured from a combi¬ 
nation of data and MC simulations. Measured efficien¬ 
cies of the high-px lepton triggers are also included in 
the overall event efficiency. Upper limits in the pro¬ 
duction cross section are calculated as a function of the 
semileptonic B.R. and shown in Figure |21 indicated by 
curved lines. For semileptonic B.R.^100%, the upper 
95% limit is 5.2 pb for a Wbb-\ike model and 2.3 pb for 
a tt model. At lower B.R., the growing limits are trun¬ 
cated because a better limit can be derived simply from 
the bigger SecVtx tagged sample. In that sample, the 
absolute normalization on the expectation of SM con¬ 
tributions other than Wbb and tt is fixed by either MC 
determination or measured from the data. The 95% C.L. 
upper limit on SM Wbb production cross section (and 
consequently any non-SM IU65-like production) is cal¬ 
culated to be 14.7 pb. Similarly, the upper limit on tt 
production, and any non-SM tt-like production, is 9.3 
pb. These are shown as horizontal lines in Figure [21 For 
comparison, the analysis in Q derived upper limits of 26 
pb(14 pb) for a IU65-like(top-like) process, assuming the 
SM semileptonic branching fraction. Moreover, based on 
the excess of events with SLT-SecVtx tags in IU±2 or 
more jets, observed during Run I at CDF and consider¬ 
ing the differences with the present analysis, we would 
expect in our dataset an excess of ~12 events due to a 
non-SM physics process. In IU±2 or more jets events, we 
exclude at the 95% C.L. a non-SM contribution exceeding 
5.6 events. 

In summary, we have compared the observed number 














of VF+jets events with at least one jet tagged by both 
the SecVtx and the SLT algorithm, with the expectation 
based on SM events and the semileptonic B.R. in heavy 
quark fragmentation. We observe no excess of events 
above expectation and set a 95% C.L. upper limit on the 
production cross section of a hypothetical non-SM parti¬ 
cle, with high semileptonic B.R. and kinematic character¬ 
istics similar to a R meson in Wbb or tt events. The limits 
are set as a function of the semileptonic B.R, providing 
the most stringent exclusion limits for the production of 
a non-SM R-like meson produced in association with a 
W boson. 
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